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Abstract 

The vortex-excitation energy on a sphere can be obtained by using the stereographic projection. 
By applying this method, we calculate the energy needed to create a vortex on a surface with a 
constant negative curvature. It is found that the energy is a linear function of the radius of the 
vortex. In accordance with this result, the interaction energy between a pair of vortices is also 
found to change linearly with the vortex separation distance. Explicit vortex configurations are 
obtained numerically with this interaction. 
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I. INTRODUCTION 



There is growing interest in topological defects on curved surfaces. One classical example 
is the Thomson problem, which addresses the question of how to configure charges on a 
sphere with minimal energy A spherical virus cell having subcellular structures on its 
surface can be viewed as a biological counterpart of the Thomson problem. The orienta- 
tional order of liquid crystal molecules in a curved sheet provides another example of the 
interaction between defects and curvature This system can be formulated as an XY 
spin model where the curvature term enters the Hamiltonian in a very similar way to that 
of the magnetic vector potential in the theory of type-II superconductors (see, for detailed 
discussions, Ref. j^j and references therein). Owing to the obvious ubiquity of spherical 
shapes, the physics of defects on a surface with positive Gaussian curvature is relatively well 
understood jj] and a number of experiments have been performed to check the theoretical 
understanding [5]. 

We have been interested in XY-type models on a surface with negative curvature 
0|. Hyperbolic geometry on such a surface is also an important model of non- Euclidean 
geometry 8j, and in physical contexts, a negatively curved surface has been introduced as 
a conceptual tool to understand disordered systems without intrinsic randomness 
To our knowledge, however, it is not entirely clear how the interaction between topological 
defects, or vortices, depends on distance in the case of negatively curved surfaces. For 



example, the potential was predicted to be short ranged in Ref. [ll| while the generalized 
Gauss law predicts it to be very long ranged Roughly speaking, the main difference 
between these two alternatives can be traced to whether or not the curvature appears as a 
source term in Gauss's law: 



E • dS = / a dV, (1) 

dV JV 

where V is a volume enclosed by a boundary surface dV . Gauss's law states that the 
surface integral of field E over dV should match with the volume integral of source terms 
distributed with density a. An important fact is that a circle on a negatively curved surface 
has an exponentially growing boundary as the radius increases. As a consequence, if only 
a fixed amount of defects contributes to the source term, the field strength should decay 
exponentially. If the curvature also serves as a source term, however, Eq. ([1]) describes 
competition between the surface and volume, both of which increase in the same exponential 



FIG. 1. Stereographic projection of a unit sphere onto a plane. A point P on a sphere, specified 
by (9,cf>), is projected onto z = re 1 ^ with r = 2 tan |. 



manner as the length scale of V grows, so the field strength never vanishes no matter how 
far away one gets from the vortex core. Although theoretical predictions about the XY 
model on a curved surface are mainly based on the former scenario {{], [llj], a numerical 
study suggests that a finite temperature is needed to unbind vortex-antivortex pairs on a 
curved surface 6|. It implies that energy should be able to compete with entropy in creating 
vortices, which will increase logarithmically with the volume of the system, meaning that 
the vortex interaction cannot decay so fast with distance. 

In this work, we directly calculate the excitation energy using the stereographic projection 
method and verify that the latter case is the correct alternative. This work is organized as 
follows: Sec. [Til reviews the stereographic projection applied to the spherical case, which 
is intended to be a mild introduction to the basic formalism. We then proceed to the 
case of negative curvature in Sec. IHIt where we also present numerically obtained vortex 
configurations based on the potential form. We conclude this work in Sec. IIVI 

II. UNIT SPHERE 

We start with reviewing how defects interact on a unit sphere, i.e., a surface with a 
constant positive curvature. The outline of this calculation has already been explained in 



Ref. 12] and the purpose of this section is to present the general method before proceeding 
to the case of negative curvature. 

Let us define coordinates (u, v) = (9,<f>) and (u',v') = (z,z) related by the stereo- 
graphic projection. As clearly seen in Fig. [IJ the projection maps the original spheri- 
cal coordinate (9, (p) onto a complex variable z = 2 tan | and its complex conjugate 



z = 2 tan I e~ % ^ . In this spherical coordinate, a point on the sphere is expressed as 
P = (X,Y,Z) = (sin 9 cos 0, sin 9 sin 0, cos 9) and the length of a line element, ds, is given 
by ds 2 = goo d9 2 + (goi + 9io) d9d(p + gn d<fi 2 = d9 2 + sin 2 9 d(f> 2 . This defines the metric 
tensor as 

„_ 9oo Soi Wl 
\9io 9n J \0 sin 2 9 
On the other hand, we want to work with the other coordinate system (V, v') and the metric 
tensor g' in the new coordinate (u', v') is transformed to g by g = U T -g' -U with the Jacobian 
matrix 
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I du' 


cV \ 
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sec 2 f 2i tan § 

sec 2 f e -i * -2i tan | 



It is now straightforward to obtain </' using the inverse matrix 



\ cos 2 § e - ** | cos 2 f e** 
icotf e"* -£catfe* 



resulting in 




= GteO • (2) 

Note that the equality between the first and second lines follows from zz = 4tan 2 | = 
4 (sec 2 | — l) . Equation ([2} means that a line element on the projected plane will be ex- 
pressed as 



ds 2 = i^dz dz 



9oo 9oi 
9w 9n 



1 



[1 + (x 2 + y 2 )/A] 2 




(3) 



if one writes z = x + iy. Hence, we can say that the determinant of the metric tensor in the 
(x, y) plane is g = det g — [1 + (x 2 + y 2 )/4] 4 . Integrating the area over the whole complex 
plane therefore yields J dS = J J yfg dx dy = J °° 2nr (1 + r 2 /4) 2 dr = Air, which is exactly 
the surface area of the unit sphere. Also note that the inverse metric tensor is given by 

9 00 9 01 \ ,/oi 

~ 9 = 10 ii = 2(1 + ^ /4) 

g 10 g 11 V l o 



where we omit the prime (') to indicate (u' , v') for brevity. When ds = w(z)\dz\, the Gaussian 
curvature K is given by the following formula js]: 
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w 2 (z) \_dzdz 
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w 2 (x, y) 
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In w(z) 
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dx 
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dy 



1 fd_ +i d_ 
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In w(x, y) 



hiw(x, y). 



(4) 



w 2 (x,y) \dx 2 dy 2 

By substituting Eq. fl3]) here, we find K = 1 for the unit sphere as expected. 

We have so far studied the fundamental property of the surface. The next step is to place 
physical objects on it. Let us consider a vector field m(x) = cos7(x) ei(x) + shi7(x) e 2 (x) 
on the surface. In tracing out its changes, however, it should be taken into account that the 
coordinate system (e^x), e 2 (x)) itself depends on the position x. A new vector field called 
the connection, or spin connection, enters here, which is derived from the given coordinate 
system by A(x) = ei(x) • Ve2(x) 13j. On the unit sphere, the most natural coordinate 
system would be obtained by differentiating P, i.e., eg = (cos 9 cos 0, cos 9 sin 0, — sin#) and 
e </> = (— sin (ft, cos (ft, 0). We express the connection in the (z,z) coordinate as (A Z ,A Z ) = 
(eg ■ de^/dz, e 9 ■ de^/dz), and it is fairly straightforward to see that A z = e e • de^/dz = 
— cos 9 (d<f)/dz). The cosine part can be easily expressed in the new coordinate system, since 
cos 9 = 2 cos 2 | — 1 = \~ zz Jt . And it follows from z = re 1 ^ that d> = Im In z = ^ (In z — In z), 
so we find that d<p/dz = (2iz)~ l and dcj)/dz = — (2iz)~ l . In short, we obtain the connection 
as 

zz/4" 



Let us now consider the contribution to the free energy due to the curvature of the surface, 
which is usually called the Frank free energy. In the one-constant approximation, that is, if 
three elastic constants associated with splay, twist, and bend are of an equal size Jjj], the 
Frank free energy assumes the following form: 

c>7 



A, 



dxp 



A t 



(6) 



with the Frank constant Ka- On a flat surface, the connection A can be set as zero and the 



remaining part describes the usual Goldstone mode. In Eq. (jfJJ), putting 7 = 
introduces one defect at 9 = and another at 9 = n. This is argued in Ref. 



] automatically 



121 ] by pointing 



out that 

if z^O, 




if z — > oo. 

The Poincare-Brouwer theorem dictates that the sum of defects on a closed surface should 
be equal to the Euler characteristic x by 

x = h I KdS = h j [v x w(x)] ' dS ' (7) 

since V x V0(x) = m5 2 (x) for a defect with charge m The Euler characteristic is related 
to the genus g of the surface, i.e., the number of handles, by \ = 2 — 2g. A sphere therefore 
has x — 2, which is of course consistent with the two defects that we have now, and the 
vector field given above connects these defects by geodesies. Additional discussion on Eq. (jSJ) 
and e lectrostatl c S is p resented in Appendix^ Si.ce ^ = ^ ' 

the calculation reduces to 



2 J (l + r 2 /4) 2 V ' ' Ar 2 \l + r 2 /AJ 



r \1 + r 2 /A t 

Let us compute the energy inside a hemisphere, which corresponds to r < 2 tan \ = 2 (see 
Fig. [1]), assuming that a defect has a very small core radius, e< 1. The answer is 

f-^jfKi^) , *-*( h i- 1 

and thus the Frank free energy over the whole sphere is given as F = 2-kKa (in - — l) . 
Alternatively, one may carry out the integration over [e, 4/e] as 

f 4/e 1 (\ -r 2 /4\ 2 (2 
F = nK A / - dr « 2?ri^ A In - - 1 



r \ 1 + r 2 /4 / \ e 

since the boundary of the defect at 9 = n is projected onto a circle of radius 4/e on the 
complex plane. Having dealt with a highly symmetric configuration as above, we may try 
a slightly more general case. We replace the defect at z = by another one at z = Zq by 
setting 7 = — Imlnz + Imln(z — Zq) and this leads to 

dl _ A 1 + 1 + 1 fl-zz/A 



dz 2iz 2i(z — z ) 2iz \l + zz/A 

1 1 / 2zz/A 



2i(z - Zq) 2iz \1 + zz/A 



FIG. 2. Projection of a unit hyperboloid onto the Poincare disk on the complex plane. The disk is 
represented by the thick solid line at Z = 0. 
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Let us denote the distance of this defect from the origin of the complex plane as \zq 
From Fig. [U we see that the image of the defect at 9 ranges over 2 tan 



P- 



2 



2 tan | ± 



esec 2 | = p±e (1 + p 2 /4) by a simple expansion and therefore the projected defect will have 



a radius of e' 
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p 2 /4). We calculate the corresponding Frank free energy 
1 
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(8) 



and the first term is evaluated as vri^^ln^- when p ~ 0(1) (Appendix [Bl). The second 
term is obtained as ^ r 2 (l + r 2 /4)~ 2 2nr dr 2ixKa (2 In 2 — l), and the last one is 
—2^7^^ In ^7 (Appendix |C|) . Summing them up, we find that 



F = 2-kKa (\^\f\ 



P74 



2ttK a In 



e cos 



In the case where two vortices are present, their configuration with minimizing the energy 
is therefore found at p = 0, where they are located at exactly opposite sides of the sphere. 



III. UNIT PSEUDOSPHERE 



In order to deal with a surface having negative curvature, first we draw a hyperboloid by 
rotating a hyperbola around its semi-major axis. The three-dimensional shape is described 
by X 2 + Y 2 — Z 2 = — 1. After parameterizing a point on the hyperboloid by P = (X, Y, Z) = 

7 



(sinh 9 cos 0, sinh 9 sin 0, cosh 9) , one may consider a projection onto a plane, Z = 0, in 
viewing the hyperboloid from Q = (0,0, —1) (Fig. [2]). By simple algebra, we find that the 
resulting point on the plane can be written as (x,y) = (r cos 0, r sin 0) with r = tanh|. 
Therefore, we have two coordinate systems, i.e., (u, v) = (#, 0) and (u',v') = (z,z) = 
(tanh | e^,tanh| e~ 1 ^), and the unit disk on the complex plane covered by this projection 
is called the Poincare disk. We want this hyperboloid to be considered as a sort of sphere , 
so we define a dot product between two vectors I = (ii,^,-^) and J = (Ji, Ji, J3) as 
I • J = I\ J\ + J 2 J 2 — I3J3. We can then simply describe the hyperboloid by P ■ P = — 1. 
According to this dot product, the line element in terms of (9, 0) is given by 

, /op op\ , 2 „ , /op dp\ JlA j2 „ , _ ; _ l2 „ 



which defines the metric tensor as 
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Note from Eq. (Q that 9 directly represents the radial distance in this metric. By using the 
Jacobian matrix 
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one can express the metric tensor in the (u', v') coordinate as 
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{9zz) , (10) 



since zz = tanh^ | = 1 — cosh - " 1 1. If one writes z = x + iy, this new metric means that 



ds 2 = {^dz dz 



9oo 9oi \ I dz 
9io 9ii I \dz 



(dx 2 + dy 2 ) , (11) 



[1 - {x 2 + y 2 )] 2 

so we find that y/g = 4 [1 — (x 2 + y 2 )} 2 on the complex plane. By using Eq. (j4j) to 
gether with Eq. f fTT]) . one can readily confirm that K = — 1. For a circle of radius R = 



8 



tanh y on the projected plane, the area inside the circle will be J dS = J J ^fg dx dy = 
/ fl 4(l - r 2 )~ 2 27rr dr = 4vr J R 2 (l - R 2 )' 1 = 47rsinh 2 Y , while it is simply ir9* 2 on the Eu- 
clidean plane. 

By differentiating the position vector P, we obtain basic unit vectors to define a coordinate 
system on this surface, i.e., e e = (cosh 9 cos 0, cosh 9 sin 0, sinh 9) and = (— sin 0, cos</>, 0). 
Since A z = e# ■ de^/dz = — cosh# |j and cosh# = 2 cosh 2 | — 1 = jzff, the connection is 
obtained as 

1 / 1 + T /, „\ 

A z = -— - = = A g . (12) 



2iz \1 — zz / 

We again consider a vector field m(x) = cos7(x) e#(x) + sin7(x) e^(x). It is notable that 
m • m = 1 is satisfied at any x with our new dot product as well. The simplest possible 
configuration would be to set 7 = to introduce a defect at (X, Y, Z) = (0, 0, 1). The Frank 
free energy is then written as 

F = ^yJJ y/9 9 aP A a A p dx a dx p 

K A f 4 ^ 9n2 1 fl + r 2 ^ 2 

x 



2 J (1 -r 2 ) 2 v ' 4r 2 \\-r 7 
Assuming a very small defect core radius e < 1, the integration up to R = tanh y yields 

2i? 2 R\ T ( ,J' , tanh^' 

ttKa ^ + In — = ttK a 2 sinh 2 — + In ■ 



l-R 2 e J \ 2 e 

Note that the first term is proportional to the total area and the second term corresponds 



to the Coulomb potential on the pseudosphere 



3: 



9* 

F c = vri^lntanhy. (13) 

In fact, we can remove the vortex at the center by assuming 7 = — Imlnz since there is no 
such restriction as Eq. (J7|) on the total sum of defects on this surface which is not closed. In 
this case, we observe 



Ka f R A r 2 



4- ^2vrr dr = 2nK A 

1 — r^Y 



R 



2 



R 2 



+ ln(l - R 2 ) 



(9* 
sinh 2 — - 9* 

The net contribution from creating the vortex is therefore 



tanh — 

AF = F - F « nK A I In ^ + 29* J . (14) 



Since the magnitude of In tanh x becomes very small as x increases, the above expression is 
approximated as AF w 2ttKa9* at any moderate distance. Interestingly, it is a logarithmic 
function of the area occupied by the vortex as in the planar case. It is reasonable to guess 
that the same functional form as in Eq. (I14p describes the interaction potential E between 
two vortices separated by a hyperbolic distance 9*, 

E{9*) « -Jq x q 2 9\ (15) 

where J ~ 2ttKa means strength of the interaction, and q\ and q 2 mean charges of the 
two vortices, respectively. The overall sign in Eq. ffl5|) is chosen so that a vortex repels 
(attracts) another vortex with the same (different) sign. If we consider the energy of a 
system containing two vortices, it can be written as H « E(9*) + Eq, where Eq is due to 
the creation of the vortices at given positions, usually given as an integral over the whole 
system. However, E may be roughly approximated as constant for a large system since 
spins far away from the vortex pair will not be much affected by small variations in their 
separation 9*. Therefore, we expect that the dominant behavior to the total energy comes 
from the interaction term, i.e., Eq. ffl5|) . while E only adds an offset. 

We can actually insert two vortices into this system, with one at the origin and the other 
away from the origin. A subtle part is that points at infinity, i.e., at \z\ = 1, should be 
equally treated since there is no reason to distinguish them. It is plausible that each field 
line should meet a point at infinity at a right angle, as does a field line emitted from a 
defect at the origin, since the exact vortex configuration near the origin will be irrelevant at 
infinity. This boundary-condition problem is indeed equivalent to that of an electric charge 
inside a conducting cylinder and one can solve this by introducing an image charge beyond 



infinity [16] . That is, for a charge at zq with \zq\ < 1, its image charge should be located at 



z 'q — Zq~ 1 in order to make every field line equally perpendicular to the boundary at infinity. 
So we need 7 = lining — zq) — Imln(z — z' ) to insert the second vortex at z = Zq, and this 
leads to 



— ~A Z 

dz 



2 



1 1 1 fl + zz 



(16) 



2i(z — Zq) 2i(z — z' ) 2iz\l — zz j 
Integrating this over a disk of radius R < 1, we indeed find that F ~ —2hKa9* + Co, where 

9* is the hyperbolic distance between the vortices and C is a system-dependent parameter 

(Appendix O). It agrees with the functional form in Eq. (fl5l) . The sign in front of 9* is 

negative due to the repulsive interaction since we are concerned with an analogy to the 
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FIG. 3. Configurations of defects on a surface of a constant negative curvature with the periodic 
boundary condition, represented on the Poincare disk, (a) When there exist two repulsive defects 
for smetic-C order, they are located at the maximum distance, (b) For hexatic order, there are 12 
defects interacting via the potential given by Eq. (|15p . (c) This pattern forms when the Coulomb 
potential [Eq. (|13p] solely comes into play, (d) Every defect has seven nearest neighbors when the 
potential is short ranged with a characteristic hyperbolic distance. 



vortex lattice in the theory of type-II superconductors. To sum up, the interaction energy 
between two vortices is a linear function of the hyperbolic distance between them unless 
they are very close. 

Before concluding this work, we briefly consider a variant of the Thomson problem, 
i.e., the Thomson problem on a pseudosphere. We can construct a periodic boundary as 



suggested in Ref. 17| in such a way that we merge every pair of opposite sides of the octagon 
in Fig. [3j The resulting closed surface has genus g = 2, and therefore \ = 2 — 2g = —2. 
The distance from one defect to another is determined by the shortest one among all the 
periodic images. Since this octagon is surrounded by 48 identical octagons, this means 
that we generally have 49 possible cases to check for determining the distance. Once the 
distance between every pair is found, it is straightforward to compute the total energy by 
using a predefined potential function and to find energy minima by applying the Metropolis 
algorithm. In terms of liquid crystals, in the smetic-C phase, the molecules are tilted when 

11 



measured relative to the surface normal and therefore described by a usual vector field. Since 
X = —2, according to Eq. ([7]), the surface should have two defects, each of which has charge 
— 1, or one dipole with charge —2. The former case is depicted in Fig. El^a), where we find 
that the two repulsive defects are located at the largest possible distance. The hexatic phase, 
on the other hand, contains sixfold orientational order and the surface has 12 defects with 
charge —1/6 each 12[. This situation is shown in Fig. [3](b) , where the defects are observed 
as pairs. Up to this point, we have used the potential given in Eq. f fT5|) . But we may also 
ask ourselves what happens with the potential shape is altered. For example, if only the 
Coulomb interaction [Eq. ( jl~3l) ] is present, a different pattern appears as shown in Fig. [3^c) . 
In addition, Fig. El^d) shows another case where the interaction potential is assumed to have 
a certain 'penetration depth' A, i.e., roughly given as exp(— 0*/X) with hyperbolic distance 
0* between a pair of defects. In Fig. [3](d) , the defects form a regular structure where each 
of the 12 defects has a coordination number k = 7 according to the Euler-Poincare relation 
and the Gauss-Bonnet the„ rem (see Ref. for details). 



IV. SUMMARY 

In summary, we calculated the vortex-excitation energy on curved surfaces by means of 
the stereographic projection. It was shown that the interaction energy is a linear function 
of the distance between vortices in the case of a negatively curved surface. This confirmed 
that the curvature should appear in the source term of the generalized Gauss law. We also 
explicitly obtained defect configurations minimizing the energy as a variant of the Thomson 
problem. 

ACKNOWLEDGMENTS 

The author is grateful to Jeong-Man Park, Petter Minnhagen, and Beom Jun Kim for 
encouraging comments. 



[1] J. J. Thomson, Phil. Mag. 7, 237 (1904). 

[2] J.-M. Park and T. C. Lubensky, Phys. Rev. E 53, 2648 (1996). 

12 



[3] A. M. Turner, V. Vitelli, and D. R. Nelson, Rev. Mod. Phys. 82, 1301 (2010). 

[4] H. Shin, M. J. Bowick, and X. Xing, Phys. Rev. Lett. 101, 037802 (2008); X. Xing, ibid. 

101, 147801 (2008); W.-Y. Zhang, Y. Jiang, and J. Z. Y. Chen, ibid. 108, 057801 (2012). 
[5] A. Fernandez-Nieves, V. Vitelli, A. S. Utada, D. R. Link, M. Marquez, D. R. Nelson, and D. A. 
Weitz, Phys. Rev. Lett. 99, 157801 (2007); T. Lopez-Leon, A. Fernandez-Nieves, M. Nobili, 
and C. Blanc, ibid. 106, 247802 (2011); T. Lopez-Leon, V. Koning, S. Devaiah, V. Vitelli, 
and A. Fernandez-Nieves, Nature Phys. 7, 391 (2011); X. Xing, H. Shin, M. J. Bowick, Z. Yao, 
L. Jia, and M.-H. Li, Proc. Natl. Acad. Sci. USA 109, 5202 (2012); J. Jeong and M. W. Kim, 
Phys. Rev. Lett. 108, 207802 (2012). 
[6] S. K. Baek, P. Minnhagen, and B. J. Kim, EPL 79, 26002 (2007); S. K. Baek, H. Shima, 

and B. J. Kim, Phys. Rev. E 79, 060106(R) (2009). 
[7] S. K. Baek, P. Minnhagen, H. Shima, and B. J. Kim, Phys. Rev. E 80, 011133 (2009). 
[8] J. W. Anderson, Hyperbolic Geometry (Springer, London, 1999). 
[9] J. C. A. d'Auriac, R. Mellin, P. Chandra, and B. Dougot, J. Phys. A 34, 675 (2001). 
[10] F. Sausset, G. Tarjus, and P. Viot, Phys. Rev. Lett. 101, 155701 (2008). 
[11] L. R. A. Belo, N. M. Oliveira-Neto, W. A. Moura-Melo, A. R. Pereira, and E. Ercolessi, Phys. 

Lett. A 365, 463 (2007). 
[12] T. C. Lubensky and J. Prost, J. Phys. II France 2, 371 (1992). 
[13] R. D. Kamien, Rev. Mod. Phys. 74, 953 (2002). 

[14] P. G. de Gennes and J. Prost, The Physics of Liquid Crystals, 2nd ed. (Oxford University 

Press, Oxford, 1993). 
[15] B. Jancovici and G. Tellez, J. Stat. Phys. 91, 953 (1998). 
[16] M. A. Furman, Phys. Rev. Spec. Top.-Accel. Beams 10, 081001 (2007). 
[17] F. Sausset and G. Tarjus, J. Phys. A: Math. Theor. 40, 12873 (2007). 
[18] F. Sausset, G. Tarjus, and D. R. Nelson, Phys. Rev. E 81, 031504 (2010). 
[19] J. D. Mancill, Am. Math. Mon. 58, 232 (1951). 



13 



Appendix A: Relation to two-dimensional electrostatics 



Let us consider the Frank free energy, given by Eq. ([6]), which has the following form 
throughout this work: 

F = 2K A JJ^(^-A.){^-A,). (Al) 

For C = (|j — A z ) — A z ), the Euler-Lagrange equation with respect to 7 is written as 

- — - 9 dC 9 dC 



d, dzd{%) dzd^L) 
<9 2 7 dA- z dA z 



(A2) 



dzdz dz dz 

Substituting A z and A z for the unit sphere [Eq. ([5])] here, we see that the last two terms 
vanish and we are left with 

<9 2 7 1 / d 2 d 2 \ 

+ ^)l(x,y) = 0, (A3) 



dzdz 4 \dx 2 dy 2 

which is Laplace's equation in two dimensions. Therefore, for a given defect configuration, 
the vector field satisfying the equivalent electrostatic problem is the one that minimizes 
Eq. flAip over spin waves. The same can be shown to be true on the Poincare disk as well 
by substituting Eq. ([12} into Eq. flM]) . 

Appendix B: Integration of the first term in Eq. 

Here we evaluate 7\ = jj D \z\~ 2 dz dz, where D is defined as a disk of \z — Zq\ < R 
containing the singularity at the origin. First, we begin with Green's theorem: 

. dM dL\ , , 
Ldx + Mdy = / / — Tr)dx dy, 



T c JjD\dx dy J 

where L and M are functions of (2, y) and C is the boundary of D. We want to find a 
vector field F = (F x ,F y ,F z ) = (L,M,F Z ) such that (V x F) 2 = f - | = (i 2 + y 2 )~ l . 
Representing this in the cylindrical coordinate, we have F = F r r + F^cf) + F z z, and rewrite 



the above expression as (V x F) 2 = \ ^(rF^) - ^ = r~ 2 . Letting F r = F z = 0, this 
yields F = r _1 (lnr + c) = r~ 2 (\nr + c) (— r sin0, r cos0, 0) with a constant c. In addition, 
we have the following identity 19j: 

I (M + iL) dz= i [M dx - L dy) + i i (L dx + M dy). 
Jc Jc Jc 
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The first part indeed vanishes for L and M above, since 



(M dx- L dy) 



fdL DM 



- 



D\dx dy 
(V ■ F) dx dy = 



dx dy 



D 



which follows from V • F = iJ^(rF r ) + + = 0. Therefore, if we define a complex 
function F(z) = L + iM as a counterpart of F, it leads to M + iL = %F and 

<L Fdz = J) (L dx + M dy) = / / (V xF) z dx dy = / / \z\~ 2 dz dz. 
Jc Jc JJd JJd 

Since F = r _2 (lnr + c)(— rsin0 — ir cos 0) = r~ 2 (lnr + c)(y + ix) = r~ 2 (lnr + c) the 
equation we are going to evaluate turns out to be 

1 



Fdz 



c 



r (lnr + c) iz dz = — 



c 



IZ 



c 



In zz + c dz 



If the contour C is given as zz = R 2 by zq = 0, then the integral becomes — § c iz~ x (In R + c) dz 
27r(lni? + c) by Cauchy's integral formula. Let us exclude a small disk 5* of radius 5^1 
around the origin to remove the constant c. The integral on the area between D and 
S is therefore 27rln-j. In the case where zq is away from the origin, we assume that it 
is on the positive real axis without loss of generality. In other words, we simply have 
z = p with < p < R — 5. The contour C is now given as (z — p)(z — p) = R 2 , or 
z = (pz + R 2 — p 2 )(z — p)" 1 . By inserting this, we get 



Fdz 



c 



c 



- In z + - h\{pz + R 2 



P 



- ln(z — p) 



dz = G(p). 



It is a bit cumbersome to directly evaluate the complex logarithms. We alternatively differ- 
entiate it with p and arrive at 



dG 
dp 



z-2p 



1 



dz. 



(Bl) 



c L 2z(pz + R 2 - p 2 ) 2z(z-p)_ 
The first term has only one pole at z = while the second has two at z = and z = p, respec- 
tively. Applying Cauchy's integral formula once again, it is found that dG/dp = —2irp(R 2 — 
p 2 y l . It is straightforward now to have G(p) = — j 2np(R 2 — p 2 )~ x dp = n In (1 — p 2 /R 2 ) + c' 
with a new constant d . However, we already know G(p = 0) = 27r(lni? + c) which deter- 
mines c'. Furthermore, we exclude S as before from the integration range. The final result 
thus becomes 

R . R 

(B2) 



Ti = 27rln- + vrln(l-p 2 / J R 2 ) = 2tt In - Jl - p 2 /R 2 . 
o o 
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If p — R — 8, for example, this formula yields T\ ~ 7rln -P, which can be cross checked by 
integrating 

r 1 ^ a r 2 -i ( m a f 1 2 -i a ~ i 2R 

/ -^2(pr dr = - cos — dr — — cos y ay ~ tt m — , 
y 5 r 2 J 5 r V2i?/ Jsy 5 

if we note that the circle centered at (-R, 0) with radius R is described as r = 2i?cos</> in 
the (r, 0) coordinate. As long as R ^> p, however, one can approximate Eq. (1B2|) simply 
as 27rln-| and inserting R = 4/e and 5 = e' gives the result in the main text. What 
happens if p > Rl Then the origin goes out of the contour and z = (p 2 — R 2 )/p comes in 
instead since p — R < (p 2 — R 2 )/p < p + R. It means that Eq. ( IB lj) now yields dG/dp = 
-2vr(p 2 + i? 2 )[2p(p 2 - i? 2 )]- 1 + 2vr(2p)- 1 = -27ri? 2 [p(p 2 - R 2 )]~\ so we obtain 

/R 2 1 
P (p 2 -i? 2 ) rfp = 71111 T^V?- (B3) 

Note that the constant of integration is determined by making this function vanish at p — > oo. 



Appendix C: Integration of the last term in Eq. ([8]) 

In this appendix, we evaluate 

T 3 = / / )~ ' ~ U J', ., dz dz. 



\z 


2 -l 


[zz + zz )/2 


\z 


- z 


2 (1 + 


z 


74) 



By rewriting this using z = r e 1 ^ and zq = p, we have 

^3= II J \ lA )rdrdc t) 

J J r 2 + p 2 — 2rpcos<p \1 + r 2 /A / 

r dr f r 2 — rp cos 4> 



l + r 2 /4 J r 2 + p 2 — 2rpcosi 
We first carry out the integration over 0: 



r 2 + p 2 — 2rpcos0 J r 2 + p 2 — rp (e^ + e -1 ^) 



(r e*P - p) (-p + r) 
i [ f I 1 



p 2 — r 2 J \e l< f > — p/r e 1 ^ — r / p 
iff! 1 



p 2 -r 2 J \( - p/r C-r/p 



d(, 



1(3 



with ( = e 1 ^. Note that we have a contour integral around a unit circle centered at the 
origin. If pjr < 1, then only the first term contributes so we get 2ir/(r 2 — p 2 ) by Cauchy's 
integral formula. If r/p < 1, on the other hand, then only the second term contributes and 
we get 27r/(p 2 — r 2 ). In short, 



/ 



dxj> 2ix 



r' 2 + p' 2 — 2rp cos </> \r 2 — p 21 



Likewise, 



COS ( 



i ( e i0 + e -i<^ 



r 2 + p 2 — 2rpcos (p J r 2 + p 2 — r p (e 1 ^ + e~ 1 ^) 

le 



L I __^±£Z__v*, 



2rp / (e*^ - p/r){e^> - r/p) ' 

_A_ / c + r 1 iC 

2rpJ (( - p/r){( -r/p) 
it f r 2 + p 2 



rp \ |r 2 — p 2 



The integration over therefore yields 



r 2 — rp cos i 



r 2 _|_ p2 _ 2rpcosi 




Gathering the terms, we see that 

" 4A rdr _ B .rui , . 2/ ,m4A , 4/e 



T 3 = 2vr / r ^ ; = 8tt [ln(l + r 2 /4)] 4A « 4vrln — - 
7 l + r 2 /4 ;Jr =" e(l 



p 2 /4)' 

which is the result in the main text. 



Appendix D: Integration of Eq. ( I16|) 



Again without loss of generality, we may set Zq = p with < p < R < 1, where i? is the 
radius of the disk over which the integration should be performed. Therefore, it follows that 
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4 = p- 1 



p' > 1. The integrand [Eq. ([IE]) ] can be then written as 

2 



1 1 1 

2i(z - p) ~ 2i(z - p') + 2iz 



1 



72 + - 
4*-P 4 



1 



+ 



1 + zz 
1 — zz 

i + k 

1 - Izl 



2\ 2 



+ 



4(z-p)(z-p') 4(z-p)(z-p') 

l /T + |d 2 \ l /l + kl 2 



4z(z — p) \1 — |z 



4^-p') vi 



Az(z — p) \ \ — \z 
1 



42(2; - p') 



1 + 


z 


2 


1 - 


z 


2 



The integration of the first term has been already done in Appendix | 

dzdz 



\<R 4 \z - p\ 

where <5 <gC 1 is the radius of a small circle around Zq to be excluded from the integration. 
We have also obtained the result for the second term: 

dzdz 7T , 1 7T . 1 



-7T = — In — 

12 4 1 



since p' > 1 [Eq. f]B3j) ]. If i? ~ 1, the contribution from these two terms will be approximately 
I In j, losing the dependence on p. By representing z = re^, we see that the third term is 
a function of r only and can be integrated directly. However, the important point is that 
the result cannot have any dependence on p. 



R 1 /l+r 2 \ 2 „ , 7T / 2 



4r 2 



lixrdr 



, In r H 

V 1 

n f R 

- In h - 

2 V 5 1 



2 



R 



r=S 



-R 2 
2R 2 



1-i? 2 , 

The fourth and fifth terms vanish together, which can be shown by a direct integration as 
follows: 

1 1 



z\<R 



+ 



dzdz 



R 



r dr 



_{z-p){z-p') (z-p)(z-p')_ 

r 2 + 1 — r cos 0(p + p' 



[r 2 + 1 — r cos 0(p + p')} + r 2 sin 2 0(p — p') 



r dr 



arctan 



r sm 1 



r 2 cos — r(p + p') + cos < 
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0. 



We proceed to the sixth and seventh terms. It is these terms that are the most relevant in 
this calculation since they describe the interaction between the two defects, with one at the 
origin and the other at z = p. Together, they can be expressed as 



1 



1 



z(z-p) z(z-p) 



I _ r 5 



1 - ^cos, 

r 



2 \^ r 2 _|_ p2 _ 2rpcos< 



Recalling Appendix (UJ we notice that 



- cos< 



■In 



7r / r 2 + p 2 



r 2 + p 2 — 2rpcos < 



\r 2 — p 2 \ 



\r 2 — p 2 \ 



2n/r 2 if r > p 
if r < p. 



(Dl) 



Therefore, the remaining integration of the sixth and seventh terms over r corresponds to 
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f R ' 1 ( 1 + r 2 \ 

/ — ( I r dr = TT \—r + 2 tanh rl 

J p r 2 \1 - r 2 J 1 J; 

= 7T [-R + 2 tanh" 1 R + p - 2 tanh" 1 p] 



(D2) 



Lastly, the eighth and ninth terms vanish according to Eq. (1D1I) since r is always smaller 
than p' > 1. 

To sum up, the p dependence essentially originates from Eq. (1D2I) . By transforming p to 
the corresponding hyperbolic distance 9* = 2 tanh" 1 p, we therefore conclude that the Frank 
free energy from the integration of Eq. (Tl6l) asymptotically results in 

2 



IK j. 



\z\<R 



dz 



A, 



dzdz w -2nK A 6* + C (R, 5), 



where Cq(R,5) is a system-dependent parameter. 
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